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Abstract 

Isoscalar monopole strength function in ^^O up to Ex — 40 MeV is discussed. We found that 
the fine structures at the low energy region up to ~ 16 MeV in the experimental monopole 
strength function obtained by the ^^0{a,a') reaction can be rather satisfactorily reproduced 
within the framework of the 4a cluster model, while the gross three bump structures observed 
at the higher energy region (16 ^ E^ ^ 4:0 MeV) look likely to be approximately reconciled by 
the mean- field calculations such as RPA and QRPA. In this paper, it is emphasized that two 
different types of monopole excitations exist in ^^O; one is the monopole excitation to cluster 
states which is dominant in the lower energy part (Ex ^ 16 MeV), and the other is the monopole 
excitation of the mean-field type such as one-particle one-hole (Ipl/i) which is attributed mainly 
to the higher energy part (16 ^ Ex ^ 40 MeV). It is found that this character of the monopole 
excitations originates from the fact that the ground state of ^^O with the dominant doubly closed 
shell structure has a duality of the mean-field-type as well as a-clustering character. This dual 
nature of the ground state seems to be a common feature in light nuclei. 



PACS numbers: 23.20.-g, 21.60.Gx, 27.20.+n 



2 



I. INTRODUCTION 



Isoscalar monopole excitation in nuclei provides important information on its underly- 
ing structure. In the collective liquid drop model, the isoscalar giant monopole resonance 
(ISGMR), which has been established in medium and heavy nuclei corresponds to a 
breathing mode of the nucleus arising due to in-phase oscillations of the proton and neutron 
fluids. In heavy nuclei, the ISGMR is observed as a single peak in the a inelastic scatter- 
ing cross sections at small angles, and its excitation energy follows an empirical formula 
Er^ ~ 80^4"^/'^ MeV, which is directly related to the compressibility of nuclear matter. A lot 
of work has been done to extract experimentally the nuclear compressibility by comparing it 
with microscopic calculations, for example, using the random phase approximation (RPA). 
Recently the isoscalar monopole distributions in ^°Zr, ^^^Sn, ^^'^Sm, and ^^^Pb were measured 
with greater precision than previously {2]. The results indicated that the compressibility of 
nuclear matter is = 231 ± 5 MeV. 

It is interesting to study what happens for the ISGMR in lighter nuclei. When the nuclear 
masses decrease from medium nuclei to light nuclei, the surface-energy correction becomes 
more important and the excitation energy of ISGMR should become lower compared with 
the empirical formula, indicating a lower nuclear compressibility |3-5|. A lot of theoretical 
work has been so far devoted to the study of ISGMR in light nuclei, for example, within the 



RPA framework 



id ] and others |11|. 



The RPA calculations with the non-relativistic framework were performed in ^^O, ^"^Ca, 
90Zr, and ^ospb Q,!!^. According to the results with the Gogny force and Skyrme forces 
etc. 4], it was found that 1) the isoscalar monopole strength in ^^O spreads out over some 
energy region of 20 ^ ^ 40 MeV with its centroid energy being = 22 ~ 29 MeV, 
the value of which depends on the A^A^ interactions employed, 2) the monopole strength 
becomes more and more concentrated in a single peak as the nucleus becomes heavier, 
3) the percentage of energy-weighted sum rule carried by the resonances increases with the 
mass number of the nucleus, and 4) the nuclear compressibility becomes smaller as the 
nucleus becomes lighter. On the other hand. Ma et al. investigated the isoscalar monopole 
modes in ^^O, ^°Ca, ^^^Zr, and ^ospb using the relativistic RPA (RRPA) method (loj. They 
found that when going from heavy to lighter nuclei, the single-peak structure of ISGMR 
in ^°^Pb changes to a peak with several small humps in ^°Ca and eventually the monopole 
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FIG. 1: Experimental isoscalar monopole strength function of ^^O [1^ is shown by the histogram. 
The experimental data below Ex ~ 10 MeV are absent because of an energy cut in the experimental 
condition. The real line is the calculated result by the relativistic RPA calculation [l^ multiplied 
by 0.25 and shifted down in energy by 4.2 MeV. This figure is taken from Ref. [l^ . 



strength spreads out widely to form a couple of peaks in ^^O. Hence, it was pointed out that 
it becomes difficult, in a nucleus like ^°Ca, to define theoretically the energy and width of 
the ISGMR. 

The experimental isoscalar monopole strengths with a great precision were recently pro- 
vided in ^^C, ^^0 and ^^Mg up to — 50 MeV, using inelastic scattering of a particles, by 
the Texas group [l^. They found that the isoscalar monopole strength in light nuclei does 
not concentrate on a single peak and the monopole strength spreads out in several regions 
of energies. The histogram in Fig. [T] shows the experimental isoscalar monopole strength 
function in ^^O [3]. It is compared with the RRPA calculation by Ma et al. [l^. It was 
found that the centroid in the RRPA response function is at 25.3 MeV, which is higher 
than the experimental data {E^. = 21.13 ± 0.49 MeV). In order to match their calculation 
to the experimental centroid, the calculated strength function was shifted down in energy 
by 4.2 MeV and furthermore they normalized it to approximately 30% of the isoscalar en- 
ergy weighted sum rule (EWSR) by multiplying the RRPA curve by a factor of 0.25 12 |. 
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FIG. 2: (Color online) RPA [dashed (black) lines] and SRPA [full ( red) lines] for the isoscalar 



monopole strength distributions of O. This figure is taken from Ref. 



Then, the normalized and shifted curve and the experimental result are in moderately good 
agreement with each other with respect to the shape of the gross three-peak structure. 
However, their calculation failed to reproduce the 0"*^ states found in the low energy region 
{5 ^ Ex ^ 16 MeV), in particular, at = 6.05, 12.05 and 14.1 MeV observed in inelastic a 
scattering and electron scattering etc. 12|, ll3|. According to the ^^0(e,e') experiments 13| . 
the three states are excited rather strongly by the (e, e') reaction, and their monopole matrix 
elements are 3.55±0.21, 4.03±0.09, and 3.3±0.7 fm^, respectively, comparable to the single- 
particle monopole strength [l^. The total percentage of the energy weighted strength to 
the isoscalar monopole EWSR (energy weighted sum rule) for these three 0"*" states amounts 

nn 

to be as large as over 15 % [13|, Il4l ]. 

In the nonrelativistic calculation for ^^O jj] a significant discrepancy is also revealed as 
compared with the experimental data, in particular, in the low energy region (5 ^ E^ < 
16 MeV), although the gross structures at the higher energy region {E^ ^ 20 MeV) in the 
RPA calculations are in rather good agreement with the experimental data. This discrepancy 
in the low-energy region can also be seen in Fig. [2] obtained by the recent second random- 
phase approximation (SRPA) calculations with a Skyrme force for ^^O [l^, in which the 
coupling between Iplh and 2p2h as well as between 2p2h configurations among themselves 
are fully taken into account. In particular, their calculation fails to reproduce the monopole 
transition strength to the 02" state at E^ = 6.05 MeV observed by the ^®0(e, e') experiment. 
Thus, the monopole strengths in the lower energy region {5 ^ E^ ^ IQ MeV) are likely to 
be out of scope in the mean field theory. These results mean that the monopole strength 
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function of ^^O is not fully understood in the mean-field theory at the present stage, and 
other degree of freedoms beyond the mean field should be taken into account. 

The 0^ and levels of ^^O including its ground state, together with their monopole 
strengths, have in the past nicely been reproduced with a semimicroscopic cluster model, 
i.e. the a+^^C orthogonality condition model (OCM) [16j. The OCM is an approximation of 
the resonating group method (RGM) Many successful apphcations of OCM are reported 
The a+^^C OCM calculation as well as the a+^^C generator-coordinate-method 



in Ref . 

one 



19| demonstrates that the 0^ state at = 6.05 MeV and the 0^ state at 



20|. 



12.05 MeV have a+^^C structures, where the a particle orbits around the ^^C(O^) core in 
an S wave and around the ^^C(2+) core in a D wave, respectively. The 14.1-MeV 0"*~ state, 

nn 

however, could not be explained by the a+^'^C model calculations |16l. Il9l|. 

Recently the structure study of ^^O has made a great advance up to E^ — 16 MeV around 
the 4a disintegration threshold. The six lowest 0"*" states of ^^O, up to E^ — 16 MeV, incluc- 
ing the ground state, have for first time been reproduced very well with the 4a OCM 
The 4a OCM shares 68 %of the energy weighted sum rule value of the isoscalar monopole 
transition of ^^O, while the a+^^C OCM does 31 %, as will be discussed below. Thus, it 
is interesting to investigate whether the 4a OCM can reproduce the experimental isoscalar 
monopole strength function in the low energy region up to E^ — 16 MeV in ^^O, a region 
which is difficult to be treated in the mean-field theory. As will be discussed below, the five 
excited 0^ states of ^^O up to E'a; ~ 16 MeV have a-cluster structures 16, 18|-20|. 

The purpose of the present paper is two fold: first is to show that the isoscalar monopole 
strength function calculated with the 4a OCM is in good correspondence to the experimental 
one in the low energy region up to E^ ^ 16 MeV shown in Fig. [H and the second is to 
emphasize two features in the isoscalar monopole excitation of ^^O, i.e. that the monopole 
excitation to cluster states is dominant in the lower energy part {E^ ^16 MeV) of the 
monopole strength function, whereas the monopole excitation of the Ipl/i-type contributes 
to the higher energy region (16 ^ -E^; ^ 40 MeV). We will show that the two features arise 
from the fact that the ground state of ^^O originally possesses a dual nature allowing a-type 
excitations as well as Iplh-tjpe ones, as will be discussed below. In this paper, a shell model 
calculation with the model space of Os-, Op-, Odls-, and 0/lp-shells for ^^O is also performed 
to investigate the extent to which the shell model works for describing the low-lying O"*" 
states. 
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In Sec. ini the monopole excitation function with the 4a OCM is formulated after a brief 
explanation of the 4a OCM framework together with the shell-model framework for ^^O. 
Results and discussions are devoted to Sec. IIIH together with the energy weighted sum rule 
of the isoscalar monopole transition. Finally we present a summary in Sec. IIVI 



II. FORMULATION 

First we formulate the isoscalar monopole strength function within the framework of the 
4a OCM. Then, the formulation of the shell model analysis is presented for ^^O within the 
model space of Os, Op, Odls, and 0/lp shells. 



A. Monopole strength function 

The strength function S{E) of the monopole excitation from the ^^O ground state 0^ is 
defined with use of the isoscalar monopole operator O = Y^ltii^i — Rcm)^ as follows, 

16 

s{E) = - ^n)l(o.tl E(^^ - ^c„.)^|o+)r (1) 

n 1=1 

where r, (i = 1 ~ 16) are the coordinates of nucleons, -Rem = jq J^lti fi is the c.o.m. co- 
ordinate of ^^O, and En denotes the excitation energy of the 0^ state of ^^O. On the other 
hand, the response function for the transition operator O is defined as 

with e representing an infinitesimal positive number. Then, R{E) is related to S{E) through 
S{E) = ~-Im [R{E)] = E ^(^ - En)mO\0t)f- (3) 

When the state |0+) is a resonance state with the complex energy En — ^r„/2, the strength 
function is expressed as 

S(E) ^ -ix,„ IRiE)] ^ i E ^^-^E^lll^lMK - Or)P, (4) 

16 

M{0^-Ot) = {0^\Y.ir^-Rcra?\Ot), (5) 

i=l 
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where r„ represents the width of the state. The isoscalar monopole transition matrix 
element, Ai{0^ — Of), has a relation with the EO transition matrix element M(EO, 0^ — Of) 
for the Oi~ and 0+ states with the total isospin T = as follows, 

16 1 _|_ 

M(EO,0+-0+) = {Oi\Y.^^ir^-Rcm)'\Of) 

1=1 

= lM{0t-0f). (6) 



The energy weighted sum rule (EWSR) of the isoscalar monopole transition [l^ reads 

J2iE^-E,)Mi0t-0f) =^xl6xi?2, (7) 

n ^ 
1 

R' = T^{Of\T.i^.-Rcm)'\Of), (8) 
1=1 

where R and m represent the r.m.s radius of the ground state and nucleon mass, respectively. 
Here, we assume that the A^A^ interaction has no velocity dependence. Employing the 



experimental charge radius of ^^O (i?c = 2.70 fm 13|), the value of R in Eq. is estimated 
to be 2.58 fm, in which the effects of the charge radius of proton ((^^)pi.oton ~ 0.8791^ fm^) 
and that of neutron (('"^)neutron ~ —0.1149 fm^) |l^ are subtracted from the charge radius 
of (R,): R = ,lRF^V)^^^^^~^V)ZZ'n = 2-58 fm. Then, the total EWSR value, 
^ X 16 X is 8.83 X 10^ fm^-MeV. 

It is instructive to see a characteristic feature of the isoscalar monopole operator in 
Eq. ([5]), which can be decomposed into two parts, internal parts and relative parts, with 
respect to 4 a clusters in ^^O (as well as a and ^^C clusters in ^^O). Since the operator 
in Eq. ([5]) has a quadratic form with respect to the coordinates of nucleons, the following 
interesting identities are realized, 

16 4 4 4 

$:(r, - R^^f = Y: E(^^+4(fc-l) - Raj' + E 4(i^a, - Ren.)', (9) 
i=l k=l 1=1 k=l 

4 16 

= $:(r, - R^)' + J2{r^ - Rc)' + 3^3^, (10) 

i=l i=5 

where Ra^ = (1/4) X^iLi ^"1+4(^-1) is the c.o.m. coordinate of the k-th a cluster, and 
(j = 1 ~ 3) are Jacobi coordinates with respect to the c.o.m. coordinates of 4 a clusters 
{Ra„ k = I ^ A): = R2- Ri, ^2 = i?3 - {Ri + i?2)/2, and = R^ - {R^ + 
R2 + Rs)/^. In Eq. ([10]), R^ = (l/4)Eti^i and Rq = (1/12) Ei=5 stand for the 
c.o.m. coordinates of a and ^^C clusters, respectively. Here we should recall the useful 
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identity of J2k=i^iRak - Rcm)^ = Y?j=il^kij^ in Eq. (i]), where /i^ (/ii = 2, /i2 = 8/3, and 
= 3) correspond to the reduced masses with respect to the Jacobi coordinates 
Equation (jO]) shows that the monopole operator consists of two parts: 1) internal parts 
[first term in the right hand of Eq. ([9])] composed of the internal coordinates of each a-cluster, 
and 2) relative parts [second term in the right hand of Eq. (j9])] acting on the relative motions 
of the 4a clusters with respect to the c.o.m. of ^^O. On the other hand, Equation f fTOj) 
shows that the monopole operator can also be decomposed into two other parts: 1) two 
internal parts, i.e. first and second terms in the right hand of Eq. ( fTOj) . composed of the 
internal coordinates of the a and ^^C clusters, respectively, and 2) relative part acting on 
the relative motion between the a and ^^C clusters. The fact that the isoscalar monopole 
operator consists of the two parts, the internal part and the relative part, plays an important 
role in the monopole excitation of ^^O, see Sec. Illli 



B. 4q OCM 

The total wave function \^(J'^) of the 4a system with total angular momentum J'^ in the 
OCM framework is expressed by the product of the internal wave functions of a clusters 
0(a) and the relative wave function \E'(J'^) among the 4a clusters 

^{r) = *(J-)0(ai)0(a2)0(a3)0(a4). (11) 

The relative wave function '^{J'^) is expanded in terms of Gaussian basis functions as 
follows, 

^{r) = Y.A.{u)^.{u), (12) 

{uF\^>{r)) = 0, (14) 

where ^2 and ^3 are the Jacobi coordinates describing internal motions of the 4a system. 
S stands for the symmetrization operator acting on all a particles obeying Bose statistics. 
u denotes the set of size parameters i^i,z/2 and 1/3 of the normalized Gaussian function, 
V^ii^y^i) = ^i.Vi'C' 6xp (— z/j,^^)y;m(05 ^ud c the Set of relative orbital angular momentum 



channels [[/i, l2]h2^ h]j depending on either of the coordinate type oi K oi H [20|, l22|, where 



Zi, I2 and Z3 are the orbital angular momenta with respect to the corresponding Jacobi 



coordinates. Equation (fT^ represents the orthogonality condition that the total wave func- 
tion f|T2l) should be orthogonal to the Pauli-forbidden states of the 4a system, up^s, which 
are constructed from Pauli forbidden states between two a-particles in OS*, OD and IS* 



states 



231]. The ground state with the dominant shell- model-like configuration (Os)'^(Op) 
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can be described properly in the present 4a OCM framework, as discussed below. 
The 4a Hamiltonian for \E'(J'^) is given as follows: 

n = j2T.-T,r. + j:[vS\^,J) + vS\^,J)_ 

i i<j 

+ E V^3a(^,J,fc)+l^4a(l,2,3,4), (15) 

i<j<k 

where Tj, V2a\i,j), V2a\'i,j), V^aihj, k) and V4q,(1, 2, 3, 4) stand for the operators of kinetic 
energy for the z-th a particle, two-body. Coulomb, three-body and four-body forces between 
a particles, respectively. The center-of-mass kinetic energy Tcm is subtracted from the 
Hamiltonian. The effective a-a interaction V2a^ is constructed by the folding procedure 
from an effective two-nucleon force. Here we take the Modified Hasegawa-Nagata (MHN) 



force 



24| as the effective NN force, which is constructed based on the G-matrix theory. It is 



noted that the folded a-a potential reproduces the a-a scattering phase shifts and energies 
of the ^Be ground state and of the Hoyle state. The three-body force is phenomenologically 
introduced so as to fit the ground state energy of the ^^C with the framework of the 3a OCM. 
The same force parameter set as used in Ref. 25| is adopted in the present calculation. In 
addition, the phenomenological four-body force is adjusted to the ground state energy of 
^^O. The three-body and four-body forces are short-range, and, hence, they only act in 
compact configurations. The coefficients Ac(z^) in Eq. (1T2|) are determined according to the 
Rayleigh-Ritz variational principle. 

The isoscalar monopole matrix element is evaluated as follows: 

16 

_^ocM(0+ _ ^ (^(0+)E(r. - R^^'mot)), (16) 

i=l 

= (^(0+)|E4(i?a, -i?cm)'|^(0|)) + 16xi?(a)2(5„i, (17) 
fc=i 

= motmi' + 1^2' + 3^/1^(0^)) + 16 X i?(a)2<5„i, (18) 
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where \1/ is the 4a OCM wave function in Eq. (1121) . In Eqs. (fT6|) ~ (fT8|) we used the relation, 



4 4 



{(f){ai)(j){a2)(j){a3)(f){a4)\^^{ri+^k^i) - ii^J V("i)</'("2)</>(a3)0(a4)) = 16 x i?(a)^j(19) 



k=l i=l 



where R{a) is the r.m.s. radius of a particle, R{a) = y\{J2i=iifi ~ RaY)a 20|, l25[. It is 
important to study the EWSR of the isoscalar monopole transition within the framework of 
the 4q! OCM. We call it the OCM-EWSR, and its definition reads 



E(i?n - E,) >lOC^(0+ - 0+) 



2 



= ^(M/(0+)|[Oocm,[H,Oocm]|^(0^)), (20) 

= — (^(0^)lE4(i?.. -i«c„.)^|vl^(Oi-)), (21) 

^ k=i 

= X 16 X (i?2 - RiaY) (22) 

Tfl/ 

where "H is given in Eq. f lT5|) and Oqcm = Y.k=i 4(-Rafe — -Rcm)^ [see Eq. ( !T7|) ]. and i? denotes 
the r.m.s. radius of the ^^O ground state given in Eq. ([8]). It is noted that the 4a OCM can 
describe the shell-model-like structure of the ^^O ground state, as shown later. Then, the 
ratio of the OCM-EWSR to the total EWSR in Eq. ([7]) is 

OCM-EWSR r^-47V_nfi8 (oo. 

total EWSR - ^ " ) -^-[2l8) " ^'^^^ 

Here we use R{a) = 1.47 fm and R = 2.58 fm, which are estimated from the experimental 



charge radii (1.68 fm and 2.70 fm, respectively 13|) with subtracting the effects of the charge 



radius of proton and that of neutron from them, the method of which is the same as that 
shown in previous section. This result means that the 4a OCM framework shares about 
70 % of the total EWSR value (the OCM-EWSR is also discussed in Appendix). This is one 
of the important reasons that the 4a OCM works rather well in reproducing the isoscalar 
monopole transitions in the low-energy region of ^^O as shown later. 

In the present paper, the energies En and isoscalar monopole matrix elements Ai in 
Eq. (jl]) are obtained by the 4a OCM calculation. As for the widths r„, we estimate the 



a-decay widths with the i?-matrix theory 29 1 



Tl = 2Pi(a)7i(a), (24) 
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7i(a)=^i(a)7^(a), (26) 



2 



-,Ua) - (27) 



where k, a and fi are the wave number of the a-^^C relative motion, the channel radius, and 
the reduced mass, respectively, and F^, Gl, and PL{a) are the regular and irregular Coulomb 
wave functions and the corresponding penetration factor, respectively. The reduced width of 
01(a) is related to the reduced width amplitude or overlap amplitude as, 01(a) = yJ^Ka), 
and the definition of yi is presented as 



4! 



3!1! 



6(r' — r) 



ydr) = -^^Y,(f')M''C) ^ m))- (28) 



Here, $^(^^0) is the wave function of ^^C, given by the 3a OCM calculation 25||, and r is 
the relative distance between the center-of-mass of ^^C and the a particle. The spectroscopic 
factor of the a+^^C(L) channel Si in the 0,^ state of ^^O, defined as 

Sl= / dr[ryL(r)]\ (29) 



is useful to analyze the obtained wave functions. 



In the present study, we perform more careful analyses than the previous ones [20|, 
in particular, for identifying the 0"*" states around the 4a threshold. The calculation of 
the resonant state in the bound state approximation is usually done by diagonalizing the 
Hamiltonian with use of a finite number of square-integrable basis wave functions. The 
positive-energy eigenstates obtained by the diagonalization are divided into resonant states 



and continuum states, and many methods for carrying out the division are proposed 30 |. 
In the present study, a pseudopotential method is adopted to divide the resonant states and 
continuum states, as shown below. 

Let's us first consider a repulsive pseudopotential V that is added to the original Hamil- 
tonian H, yielding 

H'(S) = H + SxV, (30) 

where 5 is a constant used to vary the strength of the pseudopotential. As increasing into 
negative values the constant 6 from the physical value, 6 = 0, the eigenenergy of this new 
Hamiltonian H'(6) decreases for any resonance state, which is eventually transformed into a 
bound state. On the contrary, continuum states show almost no change in their eigenvalues 

12 



as S increases into the negative region. In the present Aa OCM framework, it is important 
to study the eigenenergies with changing the constant 6 but with no change in the threshold 
energies of the a+^^C and 4a decay channels, even though we introduce the pseudopotential 
V. Here, we take the four-body potential V4Q, in Eq. f[T51) as the pseudopotential V, because 
the choice is convenient for practical reasons in the present numerical calculation. This 
pseudopotential method is simple but helpful to identify the resonant states under the bound 
state approximation. As a result, we obtained almost the same results as the previous 



ones 



20|, as will be shown below. 



Shell model calculation 



The shell model Hamiltonian of ^^O adopted here is presented as follows: 

16 16 



(31) 



1=1 i<j=l 

where ti denotes the kinetic energy of the z-th nucleon, and w*^*"-* (t;*^^^)) represents the central 
(LS) force of the effective NN interaction. The c.o.m. kinetic energy Tcm is subtracted from 
the total kinetic energy. The model space adopted here covers all configurations of Iplh and 
2p2h within the Os, Op, Odls, and Oflp shells. The spurious states of the c.o.m. motion are 
eliminated with the Lawson's method 311 ]. 

In the present study, we take the Volkov No. 2 force 32] and G3RS force 33| for v^'^^ and 
y(LS)^ respectively. The Majorana parameter (M) in the Volkov No. 2 force, the multiplying 
factor (Co^) of the G3RS, and the nucleon size parameter (6) are chosen so as to reproduce 
the total binding energies of the ground states of ^^O and ^^O, the LS-splitting between 
3/2j" and 1/2^ in ^^O, and the r.m.s. radius of the ground state of ^^O as well as possible. 
The following two parameter sets are adopted: case A for (M, (q^, b) = (0.665, 1.7, 1.70 fm) 
and case B for (0.620,2.6,2.0 fm). The present shell-model code is based on the code in 



Refs. 



34,^ 
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FIG. 3: (Color online) Comparison of ener gy spectra among experiment, the 4a OCM calcula- 



tion 



20|, and the a+^^C model calculation [16|, where the a+^^C and 4a thresholds are shown. 



Experimental data are taken from Ref. |l3l | and from Ref. 



361 ] for the state. 



III. RESULTS AND DISCUSSION 



A. 4a OCM calculation 



The energy levels of 0"*" states in ^^O obtained by the present 4a OCM calculation are 
shown in Fig. [3] and Table HI One can make the one-to-one correspondence of the six 



lowest 0"*" states observed up to • 
that the a+^^C OCM cluster model 



6 MeV in the 4a OCM calculation. It is reminded 



16| can reproduce only the lowest three 0"*" states. We 



obtained almost the same results as the previous 4a OCM calculation [20[. The six 0^ states 
have the following characteristic structures 20|: 1) the ground state (Oi") has dominantly 
a doubly-closed-shell structure, 2) the state at = 6.05 MeV and the 03" state at 
Ex = 12.05 MeV have mainly a+^^C structures 37|] where the a-particle orbits around 
the ^'^C{Of) core in an S-wave and around the ^^C(2^) core in a D-wave, respectively, 
the results of which are consistent with the previous studies with a+^^C OCM 16[ and 



the a+i^c Generator Coordinate Method (GCM) (loj, 3) the Oj {E^ = 13.6 MeV) and 
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TABLE I: Excitation energies (Ex), charge r.m.s. radii (Re), EO transition matrix elements [M(EO)], 



and particle decay widths (F) of t 



a+^^C OCM model calculation 



le 0"*" states in ^^O obtained by the 4a OCM calculation and 



16[, together with the experimental data 



13|, 



36( 1 . They are given 



in the unit of MeV, fm, fm^, and MeV, respectively. The experimental monopole matrix elements 
are obtained by the ^^0(e,e') reaction jl^. p^-™- represents the percentage of the energy weight 
strength to the isoscalar monopole EWSR [see Eq. ([7])]. The finite size effects of a particle and 



^^C are taken into account in estimating Rc with the 4a OCM and a+^^C OCM (see Ref. 
details) . 
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{E,j. = 14.1 MeV) states mainly have a+^'^C{Oi) structure with higher nodal behavior 
and a+^'^C{l~) structure, respectively, where in the latter the a particle moves around the 
i2C(l~) core (corresponding to the first 1~ state at = 10.84 MeV having an intermediate 



structure between the shell- model-like structure and cluster structure [18|) in P-orbit, and 
4) the Og state at 15.1 MeV is a strong candidate of the 4a condensate, (05*)^, with the 
probability of 61 %. 

These characteristic features of the structures of the six 0"'" states can be verified from 
the analysis of the spectoscopic factors Si defined in Eq. (12^]) . The results are shown in 
Fig. lU Since the ground state has a closed shell structure with the dominant component of 
SU(3)(A,/i) = (0,0) [26|, the values of the spectroscopic factors Si for 0^^ in Fig. HJcan be 
explained by the SU(3) nature of the state. This SU(3) character was confirmed by the recent 



no-core shell model 



27| . As mentioned above, the structures of the O2 and O3 states are 
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FIG. 4: Spectroscopic factors of the a+^^C{Ll) channels [L^ = 07, l^;, 2+, 3^;, 4+, OJ-) in the 
six 0+ sates of ^^O defined in Eq. ([29]) . 



well established as having the a+^^C(0^) and a+^^C(2i") cluster structures, respectively. 
These structures of the 0^ and 0;^ states are confirmed by the 4a OCM calculation. In 
fact, one sees in Fig. Hlthat the S"^ factors for the a;+^^C(Oj'') and a+^'^C{2i) channels are 
dominant in the and O3" states, respectively. In the O3" state, however, the S*^ factors 
of the a+^'^CiAt) channel in the 4a OCM calculation is rather large as compared with the 

n 

result of the a+^^C OCM one [16|]. This is due to the following facts: 1) The calculated 
excitation energy of the ^^C(4^) state in the present 4a OCM calculation is underestimated 
by ~ 8 MeV, while it is set to the experimental value {E^ = 14.1 MeV) in the a+^^C OCM 
calculation, and 2) thus in the 4a OCM calculation a coupled-channel effect of the a+^^C(4^) 
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channel with the a+^'^C{Of , 2f) channel is reinforced and consequently, the 5*^ factor of the 
a+^^C(4^) channel becomes larger. We expect that the S"^ factor of the a+^'^C{4f) channel 
will be smaller when the excitation energy of ^^C(4^) is properly reproduced in the 4a OCM 
calculation. 

Table [H lists the EO transition matrix elements M(EO). The M(EO) value of the 0^ state 
is reproduced well, while that of the state underestimates the experimental result, and 



this trend is similar to the result of the a+^^C OCM model 16|. On the other hand, the 
and O5" states mainly have the a+^'^C{Of) structure with higher nodal behavior and 
an a+^^C(l~) structure, respectively. The EO transition matrix element of the 0^ state 
is reproduced nicely within the experimental error (see Table H]). In Table [H the largest 
r.m.s. radius is about 5 fm for the Og" state, the wave function of which has a large overlap 
amplitude with the a+^^C(02') channel (see Fig. 3 in Ref. 20|). Hence the 5*^ factor of 
the a+^^C(02^) channel is dominant in the Og" state (see Fig. Hj), whereas those in the other 
channels are much suppressed. This dominance of the S'^ factor of the a+^^C(02') channel 
is one of the evidences for the 0^ state being the 4a-condensate-state, (05)^, because the 
Hoyle state has the main configuration, (OS")^ [25|, |39|, |40|, and the overlap amplitude of 
((05)^1 (05)^) becomes large. 

As for the decay widths of the O^^ and O5" states, the results are shown in Table HI The 
calculated width of the state is ~ 600 keV, which is quite a bit larger than that found 
for the of state ~ 200 keV. Both are quantitatively consistent with the corresponding 
experimental data, 600 keV and 185 keV, respectively. On the other hand, the decay width 
of the 0^ state is very small, 140 keV, in reasonable agreement with the corresponding 
experimental value of 166 keV, indicating that this state is unusually long lived. We should 
note that our calculation consistently reproduces the ratio of the widths of the 04", 05", and 
Og' states, i.e. about 6:3:2, respectively (see Table H]). 

Comparing the energy levels of the six O"*" states with the experimental monopole response 
function of ^^O shown in Fig. [H one notices that the energy positions of the fine structures 
in the low energy region (10 ^ 16 MeV) of the experimental response function seem 

to be in good correspondence with the energy levels of Of, Of, Of, and 0^ in Fig. [3] and 
Table [B It should be noted that the peak corresponding to the 02" state at = 6.05 MeV 
is not visible in Fig. [H because of an energy cut in the experimental condition 12| . Thus, 
it is important to study the isoscalar monopole strength function within the framework of 
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FIG. 5: (Color online) Calculated isoscalar monopole strength functions of ^^O (bold line) and 
experimental data (thin line: see Fig.[T]|l^). Here we use the calculated monopole matrix elements 
and the calculated decay widths for the six 0"'' states up to — 16 MeV obtained by the 4a OCM 
calculation (see Table H]), although the experimental excitation energies for the six O"*" states are 
employed (see text). 



the 4a OCM calculation. 



B. Isoscalar monopole excitation function with the 4a OCM calculation 

Figure shows the calculated isoscalar monopole strength function of ^^O defined in 
Eq. dl]), where we use the calculated monopole matrix elements and the calculated decay 
widths for the six O"*" states up to — 16 MeV obtained by the 4a OCM calculation (see 
Table H]), also the experimental excitation energies for the six 0"*" states are employed. We 
take into account the experimental energy resolution of 50 keV 12| for the width r„ in Eq. (j4]) 
through r„ = y^r„(OCM)^ + 0.0502, ^j^gre r„(OCM) denotes the calculated decay width 
of the n-th 0+ state of ^^O given in Table [H The calculated strength function is normalized 
so as to match the calculated strength of the 12.1-MeV peak to the experimental one. We 
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can see a rather good correspondence with the experimental data. The fine structures in 
the calculated strength function, i.e. one peak at = 12.1 MeV (corresponding to the 0^ 
state), one shoulder-like peak at E;^. = 13.8 MeV (04"), two peaks at E^ = 14.1 MeV (05") and 
15.1 MeV (Og"), are well reproduced. As mentioned above, the fine structures in the energy 
region of 10 ^ E^ ^ 15 MeV as well as the sharp peak aX E^ c:^ 6 MeV (corres pon ding to 
the state) are difficult to be reproduced by any mean-field calculations {4, 6-10, 15 1, as 
far as the present authors know. The calculated values of S{E) at each O"*" states (02" ~ Og") 
in Fig. [5] are approximately proportional to the squared values of the respective calculated 
monopole matrix elements. 

It is instructive and interesting to discuss the mechanism of why the five a cluster states 
(0^, of, O4 , 05", and 0^) of ^^O are excited relatively strongly from the ground state with 
the doubly-closed-shell-like structure [l^. Their monopole matrix elements shown in Tabled 
are comparable to the single particle strength (~ 5.4 fm^ [l^) and share about 20 % of the 
total EWSR value. Since the mechanism is closely related to the property of the ground 
state of ^^O as shown below, we first demonstrate its interesting properties with the use of 
the microscopic wave function and then discuss the monopole matrix elements in the OCM 
calculation. 

The wave function of the ^^O ground state has dominantly the doubly closed shell model 
configuration (Os)^(Op)^^ with the nucleon size parameter u = Muj/2h (M: nucleon mass), 
corresponding to the SU(3) (A, fi) = (0, 0) wave function |26[. This doubly closed shell model 
wave function is mathematically equivalent to a single cluster model wave function of a+^^C 



with the total harmonic oscillator quanta Q = 12 14| . 

1 



16! 



det|(0s)^(0p) 



12| 



X 



(-Rem)] ^ 



No] 
N2] 



'12!4! 

16! 
'12!4! 

16! 



''cm ( -Rem ) 



32z/\ 3/4 

) exp(— 16z/-R, 

TT / 



■cm J ) 



(32) 
(33) 
(34) 
(35) 



where (pcm denotes the wave function of the c.o.m. motion of ^^O, and A''o,2 are the normal- 
ization constants. (f){a) and (piC^C) stand for the internal wave function of the a cluster 
with the (Os)^ configuration and the internal wave function of ^^C with the angular momen- 



19 



turn of L belonging to the SU(3) irreducible representation (A,/i) = (0,4) for the (Os)'*(Op)* 
configuration. The relative wave functions between the a and ^^C clusters in Eqs. fl33|) and 
are described by the harmonic oscillator wave function UQimi^, f^) = '^QiiO'^Lmi^) 
with the node number n = (Q — L)/2 and Q = A. One can prove Eqs. fl5^ and fl34l) with 



help of the Bayman-Bohr theorem [28[. It should be reminded that the dominance of the 
doubly closed shell structure in the ground state of ^^O is confirmed by the no-core shell 
model with realistic A^A^ forces 

Equations (133|) and ( 134|) mean that the doubly closed shell model wave function has an 
a+^^C cluster degree of freedom. Since the monopole operator shown in Eq. flUJ]) has the 
relative part (8^3^), this can populate a+^^C cluster states by exciting the relative mo- 
tion with respect to the relative coordinate (^3) between the a and ^^C clusters. Since the 
0^ and 0;^ states of ^^O have the dominant cluster structures a+^'^C{Of) and a+^'^C{2f), 
respectively, the monopole transition strengths between the ground state and the a+^^C 
cluster states (Ot and Of) are enhanced and share 3 % and 8 % of the EWSR value, re- 
spectively 1^, IigI . fisj l . The detailed explanation using the microscopic framework is given 
in Ref. [J]. 

In the 4a OCM, the monopole matrix elements from the ^^O ground state to the a+^^C 
cluster states are evaluated with use of Eqs. (fT71) and (fT8l) . The validity of using the for- 
mulae is based on the following three facts found within the above-mentioned microscopic 
framework jl4]: 1) the ground state of ^^O is of the SU(3) (A,/i) = (0,0) nature with the 
a-clustering degree of freedom, 2) the monopole matrix elements come dominantly from the 
relative part of the monopole operator referring to the a+^^C relative motion, 8^3^, 3) the 
contribution from the other parts of the monopole operator becomes significantly smaller 
for the a+^^C cluster states, and 4) the a-cluster-type ground-state correlation significantly 
enhances the monopole strength compared with the case of the ^^O ground state being the 
pure SU(3) (A, fi) = (0, 0) wave function. These are certainly the reasons why the estimation 
of the monopole matrix element using Eqs. fll7l) and 0181) in the 4a OCM calculation gives a 
reasonable reproduction for the experimental data. In fact, the ground state of ^^O obtained 
in the 4a OCM has the dominant SU(3) (A,yu) = (0,0) component with Q = 12, and only 
the relative part with respect to the a-|-^^C relative motion, 3^3^, in the monopole matrix 
element f|T8l) gives the major contribution to the monopole matrix elements M(EO; O]*" — 0^3) 
in the present 4a OCM calculation. 
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On the other hand, the reason why the 0^ state with the 4a-gas-hke character has a 
relatively large monopole strength (~ 1 fm^) can be also understood from the property 
of the ground state of ^^O. The doubly closed shell- model wave function in Eq. f l32p is 
mathematically equivalent to a single 4a cluster wave function with Q = 12, according to 
the Bayman-Bohr theorem, as follows, 

V 16! 



. /4!4!4!4! / r r 8 

= A\ Im4o(^3>3i/) |m40(^2> 3^^)^40(^1,2//) 

X0(ai)0(a2)0(a3)0(a4)}, (36) 

where A^o is the normalization constant. This equation means that the ground state of 
^^O with the (Os)^(Op)^^ configuration inherently has a 4a-cluster degree of freedom. The 
relative part (or second term) of the monopole operator in Eq. (Q, Z]fc=i 4(-Rq,^ — -Rcm)^ = 
Yl\=i l^k$,k' 1 can excite the relative motion among the Aa particles. In other words, the 
monopole operator has an ability to populate democratically 4 a-particles by 2huj with 
respect to the c.o.m. coordinate of ^^O. The resultant state, thus, has some amount of the 
overlap with the 4a-gas-like state, i.e. Og", with the 4a-condensate-like structure [2^. The 
overlap value corresponds certainly to the monopole matrix element, M(EO). As shown in 



Ref. 



this Og state can well be described by a 4Q;-condensate-type microscopic wave 



function, called the THSR wave function 



39|]. In this THSR framework, the monopole 



matrix element to the 4a-condensate-like state is estimated to be M(EO) = 1.2 fm^, similar 
to that in the 4a OCM, M(EO; 0^^ - 0^) = 1.0 fm^ in Table [H which is calculated with 
the use of Eqs. (1171) and fllSp . Thus, the evaluation of the monopole matrix elements using 
Eqs. f[T7|) and flTH]) in the OCM framework is useful and gives a reasonable estimate for the 
monopole transition to the a-^^C cluster states and 4a-gas-like states. 

The mechanism that the 4a-gas-like state is populated by the monopole transition has a 
close connection with the mechanism that the Hoyle state with the 3a-gas-like structure is 
excited by the monopole transition, in spite of the fact that the ground state of ^^C has a 
shell-model-like compact structure with the main configuration of SU(3) (A, /i) = (0, 4). It is 



noted that the dominance o 
by the no-core shell model 



the SU(3) symmetry of the ^^C ground state 26|] was confirmed 



271 ]. This SU(3) (A,/i) = (0,4) wave function is mathematically 



equivalent to a single 3a cluster wave function, according to the Bayman-Bohr theorem, as 
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follow, 

0,=o(^^C) = |(0s)^(0p)«; (A,/i) = (0,4), r = 0+),„,_i 

0(01)0(02)0(03) I* , (37) 



~ 4!4!4! , 



«40(^2> 2^^)^40(^1,2//) 



J=0 

where iVo is the normalization factor. One can see that there exists a similarity between 
Eqs. ( I36l) and (!37j) . The latter equation means that the ground state of ^^C originally 
possesses a 3a clustering degree of freedom. As demonstrated in detail in Ref. [l^, this Set 
clustering degree of freedom in the ^^C ground state plays an important role in exciting the 
Hoyle state with the 3a-gas-like structure by the monopole transition from the shell-model- 
like ground state. The same story as in the case of the Hoyle state is also realized in the 
case of the Ojj' state with the 4a-gas-like structure in ^^0 which is excited by the monopole 
transition from the ground state with a shell-model-like structure, as discussed above. 

As for the Oc|' state, its main configuration is a+^^C(l]^) with the P-wave orbiting of 
an a cluster around the ^^C(lj~) core as mentioned above. According to the Bayman-Bohr 
theorem, the SU(3) (A, /i) = (0, 0) state of ^^O has no component of the a+^^C(l|f ) channel. 
However, the monopole strength to the 0^ state is as large as 3 fm^ (see Table H]). This is 
the reason that the state has small but important components of the a+^^C(Oi',25^) and 
0+^^0(02') configurations, as one can see from the spectroscopic factors shown in Fig.HJ It is 
noted that the a+^^C(0^) configuration is likely to be an alternative of the 4a-gas-like state 
with the dominant (OS*)^ configuration. Since these three configurations, a+^^C(05'',2^) and 
(OS*)^, can be excited from the ground state of ^^O by the monopole operator as discussed 
above, their respective contributions are coherently added to provide the relatively large 
monopole strength to the state. On the other hand, the situation of the state, 
characterized mainly by the higher nodal a+^^C(05^) state, is similar to the case of the 
0^ state. From Fig. |U the O4 state has also small but non negligible components of the 
a+^^C{Qi ,2i) and q;-|-^^C(0^) configurations, which contribute to the monopole strength 
for the O4 state. 



Shell model calculation 



Table [TTl shows the results of the shell model calculation using the spherical basis (see 
Sec. Ill Cp . The binding energies of the ground state of ^^0 and those of the l/2j" and 3/2j" 
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TABLE II: Total energies (E) and excitation energies (Ex) of the Of and states in ^^O as well 
as those of the l/2j~ and 3/2^ states in -"^^O obtained by the shell model calculation (see text). We 
also show r.m.s. radii (R) and neutron separation energies of the ^^O ground state. E, Ex, 
and Sn (R) are given in unit of MeV (fm). 
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states of are rather well reproduced in the cases A and B. We found that the excitation 
energy of the 0^ state is as large as 30 MeV in the present shell-model calculation, while 
the experimental value is as small as 6.7 MeV, the value of which is reproduced by the 
cluster model (see Table H]). It is noted that the present shell-model space includes only 
the Iplh and 2p2h configurations up to the Oflp shell. Since the Iplh configurations, 
in particular, (lsi/2)(0si/2)"\ (lp3/2)(0p3/2)~\ and (lpi/2)(0pi/2)"\ should give a crucial 
contribution to the monopole transition strengths, it is instructive to investigate how strongly 
the components of the three Iplh configurations, P{lplh), distribute among the various 0+ 
states. We found three energy regions in which the components P{lplh) become significantly 
large: 1) Ex ~ 32 MeV with P{lplh) ~ 81 %, contributed by the 0^ state, 2) Ex ~ 48 MeV 
with P{lplh) ~ 65 %, and 3) Ex ~ 60 MeV with P{lplh) ~ 40 %. This split of the Iplh 



component into the three energy region is consistent with the RRPA result [10|, ll2| in which 
the monopole strength concentrates mainly into three energy regions (see Fig. [1]). This is 
in line with the RPA and QRPA calculations (Fig. 3 in Ref. 15[), although the excitation 
energies of the three energy regions are different. The fact that the excitation energies of 
the three energy regions in the present shell-model calculations are higher than those in the 
RRPA, RPA, and QRPA calculations is reasonable, because our calculation employs the 
spherical harmonic oscillator basis, and the model space taken in its calculation covers only 
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the Iplh and 2p2h configurations within the Os, Op, Orfls, and Oflp shells (see Sec. Ill C\i . 

The fact that the present spherical shell model calculation has great difficulty to reproduce 
the low excitation energies of the Oo" and 0^ states is likely to be a common feature of the 



no-core shell model calculations 
calculations 



27(1, the FMD calculations jr-S], and the coupled-cluster 



Brown and Green 



xceptions are a few conventional shell-model works, for example, by 



46| in 1967, as far as the present authors 



42l ] in 1966 and Arima et al. 
know. Here it is instructive to briefly present their main results. 

Brown and Green discussed the low-lying three 0"^ states of ^^O with the deformed-shell 



model 



42|. It is proposed that the three 0'^ states {Of, Of, and O3") can be described by 



the mixture among the OpOh, 2p2h, and ApAh states. In their calculation, the unperturbed 
energies of the OpOh, 2p2h and 4p4h states are treated as free parameters adjusted to give the 
observed spectra, although the coupling strengths among the 0p0h,2v2h, and ipAh states 
are estimated with some approximations based on SU(3) algebra 47|. Then they found 
that the O^*" state has a dominant configuration of OpOh-tjpe, while the main configuration 
of the 0^ (of) state is of the 4p4/i-type {2p2h-tjpe). With use of the Brown-Green wave 
functions, Bertsch 43| calculated the EO transition strength between the 0^ and Of states, 
M(EO; 0^ - 05^), defined in Eq. iQ. His result is M(EO) = 1.6 - 2.8 fm^, corresponding to 
the experimental data (3.8 fm^), although the monopole transition strength between the Of 
and of states was not discussed in that paper. The important point in the Brown-Green 
calculation is that the unperturbed energy of the 4p4h configuration is taken to be lower 
than that of the 2p2h one. It is found that if the 2p2h state lies lower than the ApAh state, 
the calculated B{E2) transition rates between the resulting levels for the 0"^ states and 2~^ 
states are difficult to reconcile with the experimental data, although they could not present 
the reason why the unperturbed energy of the 4p4h configuration becomes lower than that 
of the 2p2h one. This schematic model proposed by Brown and Green was confirmed by the 
large-basis spherical shell model calculations mixing the (0-1-2-1- A)hu} excitations 4^ |45 |. 
Although they succeed in reproducing the low-lying spectrum of ^^O, the single particle 
energies are adjusted to fit six low-lying T = states in ^^O including the 02" and Of 



states 



44l |. Thus, the problem of why the excitation energy of the 0^ state is as small as 



6.05 MeV remains unclear in those shell model calculations. 

After Brown-Green's work, Arima et al. proposed a weak coupling picture 46| and showed 
that one can understand the appearance of the low- lying 0^ and Of states in ^^O if one 
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assumes weak coupling between four particles (Ap) in the sd shell and 4 hole (Ah) in p shell. 
They estimated the coupling strength with a shell model by employing the experimental 
excitation energy {Er^ = 0.86 MeV) of the l/2i state in ^^F which is described by the Ap 
in sd shell and Ih in p shell. They eventually found that the coupling strength between 
Ap and Ah in ^^O becomes significantly weak. This weakness of the Ap and Ah interactions 
is nothing but the basic assumption of a-cluster model in ^^O, i.e. one can easily imagine 
that the Ap (Ah) state corresponds to the a (^^C) cluster. The Ap in sd shell can obtain 
an extraordinarily large binding energy due to the ct-cluster correlation, and thus the states 
consisting of Ap (in sd shell) and Ah (in p shell) can be expected to lie at lower excitation 
energies, compared with the 2p2h states. 

A structure study of ^^O which explicitly treats the a-cluster degree of freedom was per- 
formed by Suzuki in 1976 with the semi- microscopic cluster model, a+^^C OCM 16j. In 
this model the relative motion between the a and ^^C clusters is solved, taking into ac- 
count the coupling between the a—^'^C relative motion and the internal rotational motion 
of ^^0(01^, 2i,Ai). Almost all levels of ^^O up to about ~ 14 MeV including the ground 
state with the dominant configuration of the doubly closed shell structure and the electro- 
magnetic transition rates (EO, E2, and E3) among them are reproduced well, together with 
the low-energy a — ^'^C scattering cross sections 48| (see Tabled for the monopole strengths). 
It is found that a lot of states in ^^O up to E^ ~ 14 MeV have a weak coupling structure 
of a+'^'^C, i.e. loosely bound a+'^'^C cluster structure. In particular the [03"] state has a 
weekly coupling structure of the a and ^"^Ci^Of) clusters [a and ^^C(2^)], in which the a clus- 
ter moves predominantly around the ^^C cluster with S-wave [Z^-wave] (see also Sec. IIII Al) . 
These successes in the a+^'^C OCM mean that 1) the weak coupling picture [46| is realized 
in ^^0, and the Ap (Ah) state of the ApAh configuration in shell model picture can be inter- 
preted as the a (^^C) cluster, and 2) the reason why the energy of the ApAh configurations 
is lower that of the 2p2h one is considered to be the a-cluster correlation for the Ap state. 
Although the a-cluster structures for the 0^ and 0^ states are much different from what 
Brown and Green thought 42|, one should mention that their conjecture, the energy of the 
ApAh configurations being lower that of the 2p2h one, is remarkable, because it has played 
an important role in the progress of our understanding the structure of ^^O. 
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D. Two features of isoscalar monopole transitions 



As discussed in Sec. IIII At the fine structures observed in the experimental monopole 
strength function in the low energy region up to Er^ ^ 16 MeV can be rather well reproduced 
within the Aa OCM framework. However, these fine structures are difficult to be reproduced 
by any mean field theory. This result means that the a clustering degree of freedom is 
inevitable to reconcile the low-energy behavior of the monopole strength function in ^^O with 
experiment. On the contrary, the RPA calculations look likely to reproduce approximately 
the three bump structures in the experimental monopole strength function in the higher 
energy region of 16 ^ -E^: ^ 40 MeV, although some normalization and energy-shifting 
procedures for the calculated strength function is needed to fit the experimental data (see 
Sec. [ID. 

Here we should remind that the ground state of ^^O is described dominantly by a doubly 
closed shell structure, (Os)^(Op)^^, in the 4a OCM calculation as well as the RPA, QRPA, 
and RRPA calculations. As discussed in Sec. IIII Al the doubly closed shell model wave 
function is mathematically equivalent to a single a-cluster wave function. This result means 
that the ground-state wave function originally has an a-clustering degree of freedom together 
with the single-particle degree of freedom. In other words, the ground-state wave function 
of ^^0 has a duality of a-clustering character and mean- field-type character. 

From these facts, one can notice that there exist two types of the isoscalar monopole 
excitation of ^^O, i.e. the monopole excitations to cluster states are dominant in the lower 
energy part {E^ ^ 16 MeV) of the monopole strength function, whereas the monopole 
excitation of the one-particle one- hole (Iplh) type contributes to the higher energy region 
(16 ^ Ex ^ 4:0 MeV). This also is in line with the first O"*" excited state of the a-particle which 
is situated at ~ 20 MeV. Thus, one can expect that the reproduction of the experimental 
isoscalar monopole strength function of ^^O in the full energy region up to E^ ~ 40 MeV will 
definitely fail, if one does not take into account simultaneously the a-cluster-type four-body 
correlations as well as the Iplh- and 2p2/i-type correlations in the structure study of ^^O. In 
order to tackle the issue, a structure calculation is desirable to be performed in which one 
uses a huge model space covering fully the a-type correlations together with the Iplh- and 
2p2h-tjpe correlations 

We here report on a trial calculation using the a+^^C cluster basis and collective basis 
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for studying the isoscalar monopole strength of ^^O 49j. In the latter respect the sym- 



plectic group Sp(6,-R) |50|, l5l| as the collective basis was used. Since the generators of the 



Sp(6,i?) group contain the monopole and quadrupole operators with respect to the nucleon 
coordinates and conjugate momenta, the group is expected to reproduce the EWSR for 
the operators. Although there is no effective A^A^ interactions which is suited for cluster- 
Sp(6,i?) mixed basis calculation, they took a phenomenological treatment, since their main 



purpose was to investigate the effect of the Sp(6,i?) group to the cluster states |52|, |53|]. They 
obtained the following two results: First comes that the effect of the Sp(6,i?) basis to the 
cluster states in the low energy region is not very large. It is noted that the and Of states 
(see Fig. [3]) are the cluster states which are reproduced by the a+^^C model. The second 
result is that there are three states in the higher energy region (20 ^ < 40 MeV) which 
correspond to the three peaks of the monopole excitation function and share about 70 % of 
the EWSR, although the excitation energies of the three states are higher by about 3 MeV. 
These results support our finding that the isoscalar monopole excitations in light nuclei, as 
mentioned above, are dominated by two features: a-cluster states at low energies and shell 
model states at higher energies. An analysis with a mixed Aa cluster and symplectic group 
basis may be useful for this study. 



IV. SUMMARY 



We have investigated the monopole strength function in the low energy region up to — 
16 MeV within the framework of 4a OCM, which has succeeded in reproducing all the six O"*" 
states observed up to E^ — 16 MeV and in giving good agreement with all of the available 
data such as the decay widths, monopole transition strengths, and r.m.s. radius of the ground 
state. It was found that the fine structures at the low energy region up to E^ — 16 MeV 
in the experimental monopole strength function obtained by the ^^0{a,a') experiment is 
rather satisfactorily reproduced within the 4a OCM framework. On the contrary, mean-field 
calculations have encountered difficulties to reproduce the fine structures of the monopole 
strength function at the low energy region as well as the monopole matrix elements for the 
0^ {E^ = 6.05 MeV), Of {E^ = 12.05 MeV), and 0^ {E^ = 14.01 MeV) states obtained 
in the ^^0(e,e') experiment. These results mean that the a clustering degree of freedom 
is inevitably necessary to reconcile the monopole strength (amounting to about 20 % of 
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the EWSR) in the low energy region with experiment. On the other hand, the 4a cluster 
model has difficulties to reproduce the gross three bump structures of the monopole strength 
function at the higher energy region of 16 ^ -E^; ^ 40 MeV. The gross bump structures look 
likely to be qualitatively reproduced by the mean-field-theory calculations such as RRPA, 
RPA, and QRPA, although the energy positions of the three bumps and the absolute values of 
the strength functions quantitatively deviate from the experimental data. In general one can 
expect the interplay between clustering degrees of freedom and mean-field degrees of freedom. 
Since the interplay affects the monopole strength function, there is a possibility that the 
isoscalar monopole strength at the higher energy region of 16 ^ ^ 40, in particular 
at its lower energy part, contains the influence by clustering degrees of freedom. The fact 
that no mean-field-theory calculations have satisfactorily reproduced the isoscalar monopole 
strength function even at the higher energy region (16 ^ ^ 40 MeV) demonstrates the 
need to investigate the interplay between clustering and mean-field degrees of freedom in 
the isoscalar monopole strength function. 

From the above results, one concludes that there exist two features of the isoscalar 
monopole excitations of ^^O, i.e. the monopole excitation to cluster states dominates the 
low energy region {E^ < 16 MeV), sharing about 20 % of the EWSR, while that to the 
Iplh-tjpe states looks likely to be predominant at the higher energy region. We indicated 
that the existence of these two types of the monopole excitations stems from the fact that 
the ground state of ^^O with the dominant doubly closed shell configuration (Os)^(Op)^^ that 
is the dominant SU(3) (A/i) = (00) symmetry has in fact a dual feature inasmuch as it can 
equivalently be described by a cluster wave function of the a-type as can be shown with the 
Bayman-Bohr theorem. When the monopole operator activates a-type degrees of freedom 
in the ground state, a-cluster states are excited, while in the case of the monopole operator 
acting on the Ipl/i-type degree of freedom in the ground state, collective states of the Iplh- 
type are populated. Thus, one will fail to reproduce the experimental isoscalar monopole 
strength function of ^^O up to E^ ~ 40 MeV if the a-cluster-type four-body correlations as 
well as the Iplh- and 2p2/i-type correlations are not simultaneously taken into account in 
the structure study of ^^O. 

The existence of two features of isoscalar monopole excitation which originates from the 
dual nature of the ground state seems to be general in light nuclei, and the case of ^^O 
discussed in the present paper is typical. This is due to the nuclear SU(3) symmetry which 
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is well verified in the ground state of light nuclei. According to the Bayman-Bohr theorem, 
an SU(3) wave function is mathematically equivalent to a cluster- model wave function. 
Thus the ground state which has a dominant SU(3) symmetry is considered to have the 
dual nature similar to the case of ^^O, which generates the two features in the isoscalar 
monopole excitation. However, this feature will be vanishing with increasing mass number, 
and, eventually only the Ipl/i-type collective motions are strongly excited, maybe, in the 
mass region beyond the /p-shell nuclei. The reason for this is that the quality of the nuclear 
SU(3) symmetry in the ground state of light nuclei is gradually disappearing because of 
the stronger effect from the spin-orbit forces in heavier nuclei. This means that the dual 
nature of the ground state is also corroding with increasing mass number. It is an intriguing 
subject to study theoretically and experimentally how these two features are changing with 
the mass number. Thus, it is strongly hoped that systematic experiments of analyzing the 
existence of these two features of monopole excitations will be performed in near future. 
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APPENDIX: Cluster Sum Rule of Isoscalar Monopole Transition 



In this Appendix, we first discuss the EWSR of the isoscalar monopole transition for 
nucleus. Then, we discuss the EWSR of the isoscalar monopole transition within the 
framework of the OCM, called the OCM-EWSR, in the case of the na OCM and two-cluster 
OCM (^^Zi and ^^Z-i with Zi = Ai/2 and Z-i = A2/2) for a self-conju gated nucleus A = An. 
Finally a general formula of the OCM-EWSR value is presented in the case of the /c-cluster 
OCM of nucleus (fc = 2, 3, ■ ■ ■), composed of the k clusters {^^Ai, ^2^2,- ■ ■, ^"Ak). 

The EWSR of the isoscalar monopole transition [l^ for nucleus [see Eq. ([7]) for ^^O] 
is given as 

A 2 



i=l 

A 



2h^ 

— xAxR\ (38) 
m 



R' = J{0t\J2i^^-Rcmy\0t), (39) 
1=1 

where R represents the r.m.s radius of the ground state, and other notations are self- 
evident. Here, we assume that the A^A^ interaction has no velocity dependence. The isoscalar 
monopole operator in Eq. fl38|) can be decomposed as 

4n n 4 n 

J2ir, - R^^f = Y: E(^^+4(fc-l) - Ra,f + E 4(i?a, - Rcraf, (40) 
1=1 k=l i=l k=l 

= E (r. - Ra.? +T.ir.- RA.f + -^^e, (41) 

where Ra^. = (1/4) X]i^=i ^j+4(fe-i) is the c.o.m. coordinate of the k-th a cluster, and Ra^ 
{RA2) stands for the c.o.m. coordinate of the "^^Zi (^^^2) cluster and £ = — Rai 
denotes the relative coordinate between the two clusters. As discussed in Sec. IllBt the 
isoscalar monopole operator in the na OCM gives non-zero contribution only for the second 
term of Eq. ( l40l) . and that in the two-cluster OCM provides with non-zero contribution only 
for the third term of Eq. fl41l) . 

In the na OCM, the total wave function of 0^ state is presented as 

n 

|0+) = $(ocM)(Q+)-Q^^^^^^ (42) 

fc=i 

where (j){ak) is the internal wave function of k-th a cluster and $('-'^'^)(0"'") stands for the 
relative wave function among the na clusters (see Sec. Ill Bp . Then, the EWSR-OCM is 
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evaluated as 

^^^(OCM) _ ^(OCM) 
P 

p 

2h^ 



4n 



{o;\j:{r.-R.n.r\ot) 

i=l 

n 



fc=l 



fc=i 

— xAx(i?2_i?(a)2)^ 

m 



(43) 



where $''°'"^^^(0p ) and E^^^^ are the eigenwave function and eigenvalue of the p-th 0"*" 
state obtained by solving the na OCM equation. Here the a-a interaction is assumed to be 
velocity- independent, and R[oi) stands for the r.m.s. radius of the a cluster. The derivation 
of Eq. (03]) should be referred to Sec. [TIE Then, the ratio of the OCM-EWSR to the total 
EWSR in Eq. ([38]) is 

OCM-EWSR _^ (R{a)~^ 
total EWSR V ^ / 

This result is common to all n and hence is the same as Eq. (l23i) for n = 4. As shown in 
Sec. [TUl the ratio for the 4a OCM is 68 %. 

In the case of the two-cluster OCM, the total wave function of O'^ state is presented as 



|0+) = $(2-clus.)^Q+)^^Ai^^)^^A2^^ 



(45) 



where (\){f^^Z\) \<\){f^'^ Z'i)\ is the internal wave function of the ^"^Z^ cluster {f^'^Z'i) and 
^(2cius.) j^g+-j g^a^j^js fQj. relative wave function between the two clusters. Then, the EWSR- 
OCM is presented as 



^^^(2clus.) _ ^( 



(2c1us.)n 



4?! 



^(^(2clus.) _ ^( 



(2clus.)^ 



9?)^ A A 

^^$(2clus.)(o+)|_li^^2|^(2clu..)(0+)), 

m Ai + A2 

2h^ 



X Ax 



m 



r2 _ ^R{A,y - ^R{A,r 



A 



(46) 



where $(2='""-)(0+) and E'J^^^^^) are the eigenwave function and eigenvalue of the p-th O"*" 
state obtained by solving the two-cluster OCM equation. Here the two-cluster potential 
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is assumed to be velocity- independent, and R{Ai) [R{A2)] stands for the r.m.s. radius of 
the ^'Zi (^^Zs) cluster. It is noted that this EWSR-OCM is realized in the case of the 
coupled- channel OCM, for example, a+^'^C{Of ,2i where there is no contribution of 

the internal monopole transitions in both the a and ^^C clusters. Then, the ratio of the 
OCM-EWSR to the total EWSR in Eq. dM} is 

0CM-EWSR(2 - cluster) _ Ai f R{Ai) y / /^(Aa) ^ 

total EWSR ~ ^^\R J R J ' ^ ' 

This ratio for the a+^^C OCM is 31 %, which is about half of that in the 4a OCM. 

Finally let us discuss the OCM-EWSR value, in general, in the case of the /c-cluster OCM 
of a ^A nucleus {k = 2,3, ■ ■ ■), composed of the k clusters {^^Ai, '^'^A^,- ■ ■, ^''Ak). The total 
wave function of 0"*" state is given as 

k 

|0+) = $(fc-cins.)(o+) J] 0(^«Z,), (48) 

i=l 

where (t){'^^Zi) is the internal wave function of the ^^Zi cluster and <|)('^~'='"'^-)(0+) stands for 
the relative wave function among the k clusters. Then, the EWSR-OCM is presented as 

An 2 

(0+|5:(r,-i?eJ'|0+) , 



^^^(fc-clus.) _ ^(fc-clus 

^ A. 



2h^ . 

— X Ax 
m 



i=i ^ 



i=l 
2 



(49) 



where R{Ai) denotes the r.m.s. radius of the ^'Zj cluster, and we assumed no contributions 
from the internal monopole transitions in the ^^Ai nucleus (z = 1,2, - ■ ■ ,k). The proof of 
Eq. (gnD is similar to those of Eqs. (gS]) and (06]). Then, the ratio of the OCM-EWSR to 
the total EWSR in Eq. ([38]) is 

OCM-EWSR(fc-cluster) _ ^ ^ Aj ( R{Ai) V 
total EWSR ~ ~ V R J ' 

The present results can be applied to the monopole transition in neutron-rich nuclei. For 
example, the ratio for the a + a + n + n OCM of ^°Be amounts to be 68 %. 
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